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ABSTRACT 

We use measurements from the South Pole Telescope (SPT) Sunyacv Zel'dovich (SZ) cluster sur- 
vey in combination with X-ray measurements to constrain cosmological parameters. We present a 
statistical method that fits for the scaling relations of the SZ and X-ray cluster observables with 
mass while jointly fitting for cosmology. The method is generalizable to multiple cluster observables, 
and self-consistently accounts for the effects of the cluster selection and uncertainties in cluster mass 
calibration on the derived cosmological constraints. We apply this method to a data set consisting 
of an SZ-selected catalog of 18 galaxy clusters at z > 0.3 from the first 178 deg^ of the 2500 deg 
SPT-SZ survey, with 14 clusters having X-ray observations from either Chandra or XMM-Newton. 
Assuming a spatially flat ACDM cosmological model, we find the SPT cluster sample constrains 
o'8(^^m/0-25)°'^° = 0.785 ± 0.037. In combination with measurements of the CMB power spectrum 
from the SPT and the seven-year WMAP data, the SPT cluster sample constrains as = 0.795 ± 0.016 
and Qm = 0.255 ±0.016, a factor of 1.5 improvement on each parameter over the CMB data alone. We 
consider several extensions beyond the ACDM model by including the following as free parameters: 
the dark energy equation of state (w), the sum of the neutrino masses (Smy), the effective number 
of relativistic species (Nca), and a primordial non-Gaussianity (/nl)- We find that adding the SPT 
cluster data significantly improves the constraints on w and Sm,y beyond those found when using 
measurements of the CMB, supernovae, baryon acoustic oscillations, and the Hubble constant. Con- 
sidering each extension independently, we best constrain w = —0.973 ± 0.063 and the sum of neutrino 
masses Sm^< 0.28 cV at 95% confidence, a factor of 1.25 and 1.4 improvement, respectively, over 
the constraints without clusters. Assuming a ACDM model with a free iVofr and Sm^, we measure 
A^eff= 3.91 ± 0.42 and constrain I]to^< 0.63 eV at 95% confidence. We also use the SPT cluster 
sample to constrain /jvz,= —192 ± 310, consistent with zero primordial non-Gaussianity. Finally, we 
discuss the current systematic limitations due to the cluster mass calibration, and future improve- 
ments for the recently completed 2500 deg^ SPT-SZ survey. The survey has detected ^500 clusters 
with a median redshift of ~ 0.5 and a median mass of ~ 2.3 x IO^^Mq/Zi and, when combined with 
an improved cluster mass calibration and existing external cosmological data sets will significantly 
improve constraints on w. 

Subject headings: galaxies: clusters: individual, cosmology: observations 
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1. INTRODUCTION 

Clusters of galaxies are the most massive collapsed ob- 
jects in the universe. Their abundance is sensitive to 
multiple cosmological parameters, in particular the mat- 
ter density, the amplitude of the matter power spectrum, 
and the dark energy equation of state (e.g., Wang & 
Steinhardt 1998; Haiman et al. 2001; Holder et al. 2001). 
Measurements of the cluster abundance that extend to 
higher redshifts become sensitive to dark energy through 
its effect on the growth of structure. This makes clus- 
ter abundance measurements an important systematic 
test of the standard dark energy paradigm, because they 
are affected by dark energy in a fundamentally differ- 
ent way than distance- redshift based tests, such as from 
type la supernovae and baryon acoustic oscillations. For 
the same reason, cluster abundance measurements also 
constrain different cosmological parameter combinations 
than distance-based tests, and their combination can 
break parameter degeneracies and achieve tighter con- 
straints than either method alone (e.g., Linder & Jenkins 
2003). 

Recently there has been significant theoretical and ex- 
perimental progress in efforts to use clusters as cosmo- 
logical probes. Large- volume numerical simulations have 
calibrated a "universal" cluster mass function over a 
broad range of cosmologies at a level better than cur- 
rent experimental uncertainties (e.g., Jenkins et al. 2001; 
Warren et al. 2006; Tinker et al. 2008; Bhattacharya et al. 
2011). Numerical simulations have also led to a better 
understanding of systematic biases in cluster mass esti- 
mates derived from a broad range of cluster observables 
(e.g., Kravtsov et al. 2006; Jeltema et al. 2008; Stanek 
et al. 2010; Becker & Kravtsov 2011). Measurements 
of the cluster abundance using optical, X-ray, and SZ 
selection methods have been used to place competitive 
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constraints on cosmology and dark energy parameters 
(e.g., Vikhhnin et al. 2d09b; Mantz et al. 2010c; Rozo 
et al. 2010; Vanderlinde et al. 2010; Sehgal et al. 2010). 
Currently, the most precise dark energy constraints from 
clusters are derived from X-ray selected samples which 
use the X-ray emission from the hot intra-cluster gas as a 
tracer of the total mass in the cluster. X-ray observables, 
particularly the gas mass and inferred pressure, tend to 
correlate with cluster mass with low scatter, indepen- 
dent of the dynamical state of the cluster or the details 
of non-gravitational physics in clusters (e.g., Kravtsov 
et al. 2006). 

Hot intra-cluster gas also causes a spectral distortion 
in the cosmic microwave background (CMB) in the direc- 
tion of clusters from inverse Compton scattering, a phe- 
nomenon known as the Sunyaev-Zel'dovich (SZ) effect 
(Sunyaev & Zel'dovich 1972). The surface brightness of 
the SZ effect is redshift-independent and largest at mm- 
wavelengths. The integrated SZ cifFect from a cluster is 
effectively measuring the cluster pressure, and is an ob- 
servable that is expected to have comparably low scatter 
with mass to the best X-ray observables (Nagai et al. 
2007; Shaw et al. 2008; Stanek et al. 2010). Therefore, a 
mm-wavelcngth SZ survey with sufficient angular resolu- 
tion is expected to provide clean, mass-limited catalogs 
out to high redshift, probing the regime where the clus- 
ter abundance is most sensitive to dark energy's effect on 
the growth rate of structure (Carlstrom et al. 2002). 

Recently, the first SZ cluster catalogs from three sur- 
veys have been released: the South Pole Telescope 
(SPT, Staniszcwski et al. 2009; Vandcrhndc et al. 2010; 
Williamson et al. 2011), the Atacama Cosmology Tele- 
scope (ACT, Marriage et al. 2011), and the Pianck satel- 
lite (Planck Collaboration et al. 2011). However, even 
with only ~10-20 clusters, the cosmological constraints 
from these surveys have been limited by the systematic 
uncertainty in the cluster mass calibration (Vanderlinde 
et al. 2010; Sehgal et al. 2011). X-ray surveys (Vikhlinin 
et al. 2009b; Mantz et al. 2010c) have achieved tighter 
constraints by adopting variations of the following cali- 
bration strategy: calibrating X-ray observable-mass rela- 
tions using X-ray hydrostatic mass estimates of relaxed 
clusters, applying this calibration to a larger sample of 
relaxed and unrelaxed clusters, and verifying the overall 
mass calibration from other methods, particularly from 
weak lensing measurements. In this work, we apply a 
similar strategy to the SPT-SZ survey using the cluster 
sample from Vanderlinde et al. (2010) (hereafter VIO), by 
incorporating an externally calibrated X-ray observable- 
mass relation and X-ray measurements of the VIO sample 
in order to present improved cosmological constraints. 

This paper is organized as follows. In Section 2, we 
describe the relevant SZ, X-ray, and optical data, anal- 
ysis methods, and the external cosmological data sets 
used in this work. In Section 3, we describe and imple- 
ment a self-consistent cosmological analysis using SZ and 
X-ray observations of the SPT cluster sample that siiinil- 
taneously constrains cosmology and the relevant SZ and 
X-ray cluster scaling relations while accounting for the 
SPT cluster selection function. In Section 4, we discuss 
the constraints on a ACDM cosmological model from the 
SPT cluster sample, and compare our results to the con- 
straints from observations of the CMB power spectrum. 
In Section 5, we consider extensions to the ACDM model 
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by including the following as free parameters: dark en- 
ergy equation of state, the sum of the neutrino masses, 
the effective number of relativistic species, and a primor- 
dial non-Gaussianity. We report the relative improve- 
ments using the SPT data to constrain each extension. 
In Sections 6 and 7, we discuss the limiting systematics 
and implications for applying this method to the larger 
SPT cluster sample. 

In this paper, unless otherwise specified, the cluster 
mass will refer to M500, the mass enclosed within a 
spherical radius, rsoo, where the cluster's mean matter 
density is 500 times the critical density of the universe 
at the observed cluster redshift. The critical density is 
Pcrit{z) = 3H^{z)/8'itG, where H{z) is the Hubble pa- 
rameter. 

2. DATA AND OBSERVATIONS 
2.1. Cluster Data and Observations 

The cluster sample used in this work is a sub-sample of 

a SZ-selected catalog from the SPT that was described 
in VIO. The VIO catalog consisted of 21 clusters selected 
by their SZ significance from 178 deg^ of sky surveyed by 
tiic SPT in 2008. As in VIO, wc use only the 18 clusters 
at z > 0.3 for the cosmological results in this work. The 
optical and X-ray properties of this catalog have been 
described previously in High et al. (2010) and Anders- 
son et al. (2011), hereafter HIO and All, respectively. 
In this section, wc summarize the VIO, HIO, and All 
data sets, analysis, and results used in this work. We 
also report additional spectroscopic redshift and X-ray 
measurements for several clusters. 

2.1.1. SZ Observations and the Cluster Sample 

The lO-mcter diameter SPT is a mm-wavelength tele- 
scope designed to conduct a large-area survey with low 
noise and ~1 arcminute angular resolution. The SPT re- 
ceiver consists of a 960 element bolometer array that is 
sensitive in three bands, at 95, 150, and 220 GHz. De- 
tails of the telescope and receiver can be found in Padin 
et al. (2008), Carlstrom et al. (2011), and Dobbs et al. 
(2011). The primary goal of the SPT survey is to search 
for clusters of galaxies via the SZ effect in a 2500 deg^ 
survey that was completed in November 2011. 

The first cosmological constraints from the SPT clus- 
ter survey were reported in VIO, with an accompanying 
cluster catalog. These results were derived from SPT 150 
GHz observations of 178 deg^ observed in 2008, from two 
approximately equal area fields centered at right ascen- 
sion (R.A.) 5''30'", declination (deck) -55° and R.A. 
23''30'", deck -55°. Cluster candidates were identified 
in the SPT maps by using a matched spatial filter tech- 
nique (Haehnelt & Tegmark 1996; Melin et al. 2006). In 
brief, the SPT maps are filtered in Fourier space to opti- 
mize the detection of cluster-like objects using a source 
template constructed from a /3-model of variable angu- 
lar size. This is done while accounting for the expected 
signals from the dominant sources of astrophysical con- 
tamination, instrumental and atmospheric noise, and the 
effects of the SPT beam and timestrcam filtering. Can- 
didate galaxy clusters were assigned an SZ significance, 
^, defined as the highest signal-to-noise across all filter 
scales. 

VIO used simulations to characterize the SZ selection 



function and the scaling between ^ and cluster mass. 
Simulated SZ maps were generated from large-volume 
dark matter simulations (Shaw et al. 2009) using the 
semi-analytic gas model of Bode et al. (2007). The 
gas model was calibrated to match the observed X-ray 
scaling relations for low-redshift {z < 0.25) clusters. 
The cluster selection was characterized by applying the 
matched filter to multiple sky realizations that included 
the dominant astrophysical components (primary and 
lensed CMB, thermal SZ, and point sources), instru- 
mental and atmospheric noise, and the SPT filtering. 
These simulations found that at ^ > 5, the SPT cata- 
log was expected to be ^ 95% pure. This result is con- 
sistent with optical follow-up which found optical clus- 
ter counterparts to 21 of the 22 candidates above this 
threshold. The 21 optically confirmed clusters had a me- 
dian redshift oi z = 0.74, and the sample was predicted 
to be nearly 100% complete above a mass threshold of 
M500 6x lO^/i^^Af© at z = 0.6. The simulations were 
also used to put conservative priors on the ^-mass rela- 
tion in the VIO cosmological analysis. Even with only 18 
clusters, the improvement in the cosmological constraints 
was limited by the assumed systematic uncertainty on 
the normalization of the ^-mass relation. 

The full cluster catalog used in this work is given in Ta- 
ble 1. For each cluster, we report the name, position, red- 
shift, and the SZ and X-ray observables, where the latter 
assumes a default cosmology. We note that the only SZ 
product needed for the cosmological analysis described in 
Section 3 is the SZ observable ^. In this work we improve 
the cosmological constraints relative to VIO by reducing 
the uncertainty on the ^-mass relation through inclusion 
of X-ray observables which have an obscrvablc-mass rela- 
tion that has been externally calibrated, as described in 
Vikhlinin et al. (2009a) and summarized in Section 3.1.1. 

2.1.2. Optical Redshifts 

Redshifts of the SPT clusters were measured through 
a combination of optical photometry and spectroscopy. 
The majority of the observations and data analyses are 
identical to those in HIO, to which we refer the reader for 
a more detailed description. Relative to HIO, we include 
spectroscopic redshift measurements for seven additional 
clusters, which we briefly describe here. All cluster red- 
shifts are given in Table 1. 

Optical counterparts and photometric redshifts were 
measured from a combination of imaging from the Blanco 
Cosmology Survey (BCS, see Ngeow et al. 2009) and tar- 
geted observations using the Magellan telescopes. Opti- 
cal images were searched for red sequence galaxies within 
a 2' radius of the SPT candidate location. A cluster was 
identified through an excess of red sequence galaxies rela- 
tive to the background, and the photometric redshift was 
estimated by fitting a red sequence model. The redshift 
uncertainty varies over the sample, however it is typically 
Az/{l + z) - 0.03. 

For 15 of the 18 clusters, wc have also measured spec- 
troscopic redshifts, which we use for the cluster's red- 
shift when measured. For eight of the clusters we use 
the spectroscopic redshifts as reported in HIO, which 
were measured using the Low Dispersion Survey Spec- 
trograph (LDSS3) on the Magellan Clay 6.5-m telescope. 
For SPT-CL J0546-5345, we use the redshift reported in 
Brodwin et al. (2010), measured using multi-slit spec- 
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TABLE 1 

The SPT 178 deg^ Cluster Catalog and Observables 



Object Name 



R.A. 

(deg) 



decl. 
(deg) 



Photo-z Spec-z 



Yx 



(10^^ Mq keV) 



SPT 


CL 


J0509-5342 


77.336 


-53 


705 





47(4) 


0.463 


6.61 


4.3 ±0.8 


SPT- 


CL 


J0511-5154 


77.920 


-51 


904 





74(5) 




5.63 




SPT- 


CL 


J0521-5104 


80.298 


-51 


081 





72(5) 




5.45 




SPT- 


CL 


J0528-5300 


82.017 


-53 


000 





75(5) 


0.765 


5.45 


1.6 ±0.5*= 


SPT- 


CL 


J0533-5005 


83.398 


-50 


092 





83(5) 


0.881 


5.59 


1.0 ±0.4'' 


SPT- 


CL 


J0539-5744 


85.000 


-57 


743 





77(5) 




5.12 




SPT- 


CL 


J0546-5345 


86.654 


-53 


761 


1 


16(6) 


1.067'' 


7.69 


4.8± 0.8*^ 


SPT- 


CL 


J0551-5709 


87.902 


-57 


156 





41(4) 


0.423 


6.13 


1.9 ± 0.4'' 


SPT- 


CL 


J0559-5249 


89.925 


-52 


826 





66(4) 


0.611 


9.28 


6.4 ±0.8 


SPT- 


CL 


J2301-5546 


345.469 


-55 


776 





78(5) 


0.748'' 


5.19 




SPT- 


CL 


J2331-5051 


352.958 


-50 


864 





55(4) 


0.571 


8.04 


3.5 ±0.6 


SPT- 


CL 


J2332-5358 


353.104 


-53 


973 





32(3) 


0.403" 


7.30 


6.1 ± 0.8'' 


SPT- 


CL 


J2337-5942 


354.354 


-59 


705 





77(5) 


0.781 


14.94 


8.5 ± 1.7 


SPT- 


CL 


J2341-5119 


355.299 


-51 


333 


1 


03(5) 


0.998 


9.65 


4.7 ± 1.0 


SPT- 


CL 


J2342-5411 


355.690 


-54 


189 


1 


08(6) 


1.074" 


6.18 


1.4 ±0.3'' 


SPT- 


CL 


,12355-5056 


358.955 


-50 


937 





35(4) 


0.320" 


5.89 


2.2 ±0.4'' 


SPT- 


CL 


J2359-5009 


359.921 


-50 


160 





76(5) 


0.774" 


6.35 


1.8 ±0.4'' 


SPT- 


CL 


JOOOO-5748 


0.250 


-57 


807 





74(5) 


0.701" 


5.48 


4.2 ± 1.6'' 



Note. — ,f is the maximum signal-to-noisc of the SPT-detection obtained over the set of filter scales for each cluster. The 
cluster positions in R.A. and decl. are given in degrees and refer to the center of the SZ brightness in the SPT map filtered 
at the preferred scale to maximize the signal-to-noise. We give the estimated photometric redshift and spectroscopic redshifts, 
where available. To be consistent with All, Yx is calculated assuming a preferred ACDM cosmology using WMAP7±BAO±-f/o 
data with Om = 0.272, = 0.728 and Ha = 70.2 km s"' Mpc"' (Komatsu at al. 2011). In §4 and §5, Yx is recalculated as a 
function of cosmology and scaling relations for each step in the Markov chain. 

"New spectroscopic redshift since VIO. 
''Updated Yx since All. 



troscopy with the Inamori Magellan Areal Camera and 
Spectrograph (IMACS) on the Magellan Baade 6.5-m 
telescope. Finally, there are six clusters that have new 
spectroscopic redshifts, which we report in this work for 
the first time in Table 1. These redshifts were measured 
with a combination of IMACS and GMOS on Gemini 
South, and the details of the data and analysis will be 
described in Bazin et al. (2012, in prep.). 

2.1.3. X-ray Observations 

X-ray observations were obtained using Chandra and 
XMM-Newton for 14 of the clusters in Table 1. The 
majority of the X-ray observations, data reduction, and 
analyses are the same as described by All, to which we 
refer the reader for a more detailed description. Rela- 
tive to All, we include new Chandra observations for 
five clusters, and re-run the X-ray analysis for the five 
clusters with new optical spectroscopic redshifts, one of 
which also had new Chandra observations. In this sec- 
tion, we summarize the X-ray observations and results, 
and describe additional analyses required to incorporate 
the X-ray measurements in the cosmological analysis. 

Summarizing All, 15 of the 16 highest ^ clusters from 
VIO were targeted for X-ray observations, however, in 
this work, we use only the 14 clusters at z > 0.3. Of 
these, twelve were observed with Chandra and four clus- 
ters were observed with XMM-Ncwton. Two clusters 
were observed by both Chandra and XMM-Newton, and 
for these clusters only the Chandra data was included 
in the analysis. From the data, the X-ray observables, 
Mg, Tx, and Yx, were measured in a manner identical 



to Vikhlinin et al. (2009a), where Mg is the gas mass 
within rsoo, Tx is the core-excised X-ray temperature in 
an anmihis between 0.15 — 1.0 x r5oo, and Yx = MgTx- 
We solved for each observable and r5oo iteratively, to 
maintain consistency with their respective observable- 
mass relations. Since All, five of the clusters have new 
spectroscopic redshift measurements. For these clusters 
we repeat the All reduction and analysis using the new 
redshifts, and give the updated resiilts in Table 7. 

Five of the chistcrs from All have had additional 
Chandra observations, which we include in this work. In 
Table 8, we list these clusters, the Chandra observation 
IDs, and the improvement in exposure time and cluster 
source counts adding the new observations. We repeat 
the All reduction and analysis to derive new constraints 
on the X-ray observables, which are given in Table 7. 
For these results, relative to All, we use more recent 
Chandra analysis software (CIAO 4.3) and calibration 
files (CALDB 4.3.3). We find that the new Chandra cal- 
ibration files typically change Yx by < 5%. This is at a 
level below the assumed mass-normalization imcertainty 
that we assign in Section 3.1.1. 

For the cosmological analysis in this work, described 
in Section 3, we need to calculate the X-ray observables 
as a function of cosmology and scaling relation param- 
eters. To do this, we derive density and temperature 
profiles for all 14 clusters with X-ray data. We calculate 
Tx{r) and Mg{r) (for the calculation of Yx{r)) from the 
X-ray observations of each cluster assuming a reference 
cosmology, where r corresponds to a physical radius in 
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the cluster and the profiles are defined to return the clus- 
ter observable within r. The reference cosmology is cho- 
sen to match All; a preferred ACDM cosmology using 
WMAP7+BAO-FiJo data with SIm = 0.272, = 0.728 
and Ho = 70.2 km s'^ Mpc-^ (Komatsu et al. 2011). 
For three of the clusters with the lowest X-ray photon 
counts, the Tx{r) profiles have jumps which appear un- 
physical. For this reason, we have assumed a functional 
form of Tx{r) that, in combination with the measured 
Mg{r), matches the pressure profile from Arnaud et al. 

(2010) and is normalized to give the measured Yx as- 
suming the reference cosmology. When considering the 
eleven clusters with well-behaved tempcratiirc profiles, 
we find that our cosmological results in Section 4 negli- 
gibly change when assuming either the functional form 
of Tx(r), or the profile derived from the data. Therefore, 
we consider this approximation valid for this work. 

2.2. External Cosmological Data Sets 

In addition to the SPT cluster data set, we incorporate 
several external cosmological data sets, including mea- 
surements of the CMB power spectrum (CMB), the Hub- 
ble constant (Hq), baryon acoustic oscillations (BAO), 
type la supernova (SNe), and big bang nucleosynthesis 
(BEN). We will use these abbreviations when referring 
to these data sets, and will use several different combi- 
nations of them in our analysis and results. Below we 
give references and a brief description of each external 
data set. Also, when discussing our results in Section 4 
and onward, we will define the SPTcl data set as the 
combination of the SPT-SZ data, optical redshift, and 
X-ray measurements described in Section 2.1. 

We use measurements of the CMB power spectrum 
from the seven-year WMAP data release (WMAP7, Lar- 
son et al. 2011) and 790 deg^ of sky observed with 
the SPT (Keisler et al. 2011). Following Keisler et al. 

(2011) ,^ we fit the CMB data to a model including pri- 
mary CMB anisotropy plus three nuisance parameters 
that model "foreground" signals detectable in the SPT 
data. We use low-redshift measurements of Hq from 
the Hubble Space Telescope (Riess et al. 2011), which 
we include as a Gaussian prior of Hq = 73.8 ± 2.4 
km s~^ Mpc~^. We use measurements of the BAO 
feature using SDSS and 2dFGRS data (Percival et al. 
2010). The BAO constraints have been applied as a 
measurement of rs/Dv{z = 0.2) = 0.1905 ± 0.0061 and 
rs/Dv{z = 0.35) = 0.1097 ± 0.0036; where r, is the 
comoving sound horizon size at the baryon drag epoch, 
Dv{z) = [{l + z)'^D\{z)cz/H{z)]^/^, Da{z) is the angu- 
lar diameter distance, and H(z) is the Hubble parameter. 
The inverse covariance matrix given in Eq. 5 of Percival 
et al. (2010) is used for the BAO measurements. We 
use measurements of the luminosity distances of Type la 
supernovae (SNe) from the Union2 compilation of 557 
SNe (AmanuUah et al. 2010), and include their treat- 
ment of systematic uncertainties. Finally, we use a BBN 
prior from measurements of the abundances of ^He and 
Deuterium (Kirkman et al. 2003), which we include as a 
Gaussian prior of Ofc/i^ = 0.022 ± 0.002. 

3. COSMOLOGICAL ANALYSIS 



In this section, wc outline the cosmological analysis 
method for the SPT cluster data set, including the cal- 
culation of the cosmological likelihood and the assumed 
parameterization for the cluster mass-observable rela- 
tions. This implementation allows for self-consistent 
constraints on cosmology and the cluster scaling rela- 
tions, i.e., the chistcr mass calibration, by simultane- 
ously varying the cluster-mass observable relations and 
cosmological parameters using a Markov Chain Monte 
Carlo (MCMC) technique. The method is generalizable 
in way that can include additional cluster observables 
from other data sets (e.g., weak lensing shear, optical ve- 
locity dispersions) . We have incorporated our calculation 
of the SPT cluster hkclihood into the CosmoMC code^ 
of Lewis & Bridle (2002) to compute its joint likelihood 
with the external cosmological data sets. 

3.1. Scaling Relation Parameterization 
3.1.1. X-ray: Yx - Msoo 

Following Vikhlinin et al. (2009b), we use Yx as an 
X-ray proxy for cluster mass, M500. We assume a 
Yx — M500 relation of the form 



500 



1014Mq//i 



= (Axh'/') 



Y: 



X 



3 X IQi^MfT, keV 



Eiz) 



Cx 



(1) 

parameterized by the normalization Ax, the slope Bx, 
the redshift evolution Cx, where E{z) = H(z)/Hq, and 
a log- normal scatter Dx on Yx- Relative to the form of 
this equation in Vikhlinin et al. (2009b), we have multi- 
plied the right-hand side by an extra factor of h, so that 
the cluster mass M500, is in units of Mq/H to match the 
( — M500 relation in Section 3.1.2. For our cosmological 
analysis, we assume Gaussian priors on the scaling re- 
lation parameters, which wc list in Table 2. The priors 
arc motivated by constraints from X-ray measurements 
by Vikhlinin et al. (2009a) and simulations, which we 
describe below. 

Vikhlinin et al. (2009a) constrained the Yx — M500 re- 
lation using X-ray observations of a low-redshift (z < 0.3) 
sample of 17 relaxed clusters to estimate the hydrostatic 
total mass and Yx- From simulations, Kravtsov et al. 
(2006) put an upper limit on the systematic offset in 
the Yx — M500 relation between relaxed and unrelax;ed 
clusters of 4%. Simulations also expect that biases in 
hydrostatic mass estimates arc less for relaxed clusters 
and typically < 15% (Nagai et al. 2007). Therefore, a 
Yx — M500 relation calibrated from hydrostatic mass es- 
timates of a relaxed cluster sample should have minimal 
biases and be appUcable to a larger cluster sample of 
both relaxed and unrelaxed clusters. 

From the above measurements, Vikhlinin et al. (2009a) 
obtained best-fit values of Ax = 5.77 ± 0.20 and Bx = 
0.57 ± 0.03, where the uncertainties are statistical only. 
They estimated the systematic uncertainty in the Ax 
calibration by comparing to weak lensing mass estimates 
from Hoekstra (2007) for a sample of 10 low-redshift clus- 
ters. From this analysis, they estimated a la uncertainty 
of 9% on the Chandra mass scale calibration, which we 
add in quadrature with their quoted statistical uncer- 
tainty on Ax ■ We have therefore assumed Gaussian pri- 



^ http://pole.uchicago.edu/public/data/keislerll 



^ http: / /cosmologist.info/cosmomc/ 
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ors of Ax 5.77 ± 0.56 and Bx = 0.57 ± 0.03. 

Wc assume a Gaussian prior of Cx = —0.4 ± 0.2, con- 
sistent with self-similar evolution and a 50% uncertainty. 
There are relatively few observational constraints on the 
normalization of the Yx — M^qq relation at ^; > 0.3. 
However, this level of Cx uncertainty was chosen to 
match the prior assumed in Vikhlinin et al. (2009b), 
which was in turn motivated by constraints from the sim- 
ulations of Kravtsov et al. (2006). The la prior on Cx 
would correspond to a 6% difference in the mass calibra- 
tion between z = 0.0 and 0.6. For the highest redshift 
cluster in this work, z = 1.074, this would correspond to 
an additional 11% uncertainty in the mass calibration, 
and a 15% total uncertainty when also considering the 
prior on Ax- 

We assume a Gaussian prior of Dx = 0.12 ± 0.08, 
which we truncate below 0.02 and where Dx = 0.12 cor- 
responds to a 12% log-normal scatter in Yx for a given 
mass. This scatter has been measured to have values 
ranging from 0-12% (Vikhlinin et al. 2009a; Mantz et al. 
2010a). Analogous to Vikhlinin et al. (2009b), for our 
cosmological analysis we have chosen a prior centered 
on a value which is consistent with simulations (e.g., 
Kravtsov et al. 2006). Although we have assumed a 
larger uncertainty on the scatter than the range typi- 
cally found in simulations, this uncertainty is negligible 
for our cosmological constraints, see Section 6. 

3.1.2. SZ:C-Mr,QQ 

As in VIO, we use the detection significance ^ as an SZ 
mass proxy. However, since the relation between ^ and 
halo mass is complicated by the comparable effects of 
intrinsic scatter and instrumental noise, we introduce the 
unbiased significance, (: the average detection signal-to- 
noise of a simulated cluster, measured across many noise 
realizations, evaluated at the preferred position and filter 
scale of that cluster as determined by fitting the cluster 
in the absence of noise. 

We relate unbiased significance C to the detection sig- 
nificance £, as follows. Firstly, C is related to (^) through 
the relation 

C=./W^ (2) 

at ^ > 2. This maximization bias comes from having 
maximized ^ across possible cluster positions and filters 

scales, effectively adding three degrees of freedom to the 
fit with ^ analogous to a x^- Additionally, {(_) relates to 
^ by a Gaussian scatter of unit width. Simulations have 
been used to verify that these approximations introduce 
negligible bias or scatter compared to the Poisson noise 
of the sample. For further details we refer the reader to 
VIO. 

We assume a C — Afsoo relation of the form 



C = A 



SZ 



500 



3 X lOi^Mn/i- 



Bsz 



E{z) 
£;(0.6) 



Cs2 



(3) 



parameterized by the normalization Asz, the slope Bsz, 
the redshift evolution Csz, and a log-normal scatter, 
Dszi on C- VIO motivated the form of this relation 
based on physical arguments, and the expected range 
of these parameters based on self-similar arguments. In 
VIO, the cluster mass was defined within a spherical re- 
gion in which the density is equal to 200 times the mean 



matter density at the cluster redshift. In this work, to be 
consistent with the Yx — -W500 relation, we are defining 
the cluster mass as M500, the mass in a spherical radius, 
Tbooi within which the density is equal to 500 times the 
critical density of the universe at the cluster redshift. 
This change has motivated a change in the redshift evo- 
lution term from (1 + z) to E{z), because of the expected 
self-similar scaling between Ysz and M500 (e.g., Kravtsov 
et al. 2006) . In addition, we allow for a correlated scatter 
between ( and Yx with a correlation coefficient p, which 
we allow to uniformly vary between 0.02 and 0.98, but 
away from and 1 for numerical reasons. 

Analogous to VIO and summarized in Section 2.1.1, 
we have used simulated SZ maps to characterize the 
scaling between ( and cluster mass. We have repeated 
this exercise to match the form of the scaling given in 
equation 3, and we give the Gaussian priors in Table 2. 
The fractional uncertainty on each parameter matches 
VIO, except for the log-normal scatter, for which we al- 
low a larger uncertainty in this work. However, this un- 
certainty remains negligible for these cosmological con- 
straints, see Section 6. 

3.2. Likelihood Model 

The analysis method employed in this work closely mir- 
rors the one presented by VIO with extensions to incor- 
porate the X-ray data. In VIO, the parameter space was 
explored through importance sampling of pre-existing 
WMAP MCMC chains. In this work, we have elected 
to utilize a full MCMC algorithm. This is accomplished 
through the use of the CosmoMC analysis package, where 
we have included the cluster abundance likelihood as an 
additional module in the CosmoMC likelihood calcula- 
tion. Among the numerous advantages to this approach 
is the ability to enforce quantitative convergence crite- 
ria as well as the optional inclusion of supplemental data 
sets. 

Each step in the Markov chain selects a new point 
in the joint cosmological and scaling relation parame- 
ter space. Prior to passing these variables to the cluster 
likelihood evaluation, we use the Code for Anisotropies in 
the Microwave Background (CAMB) (Lewis et al. 2000) 
to compute the matter power spectrum at 20 logarithmi- 
cally spaced redshifts between < z < 2.5. The matter 
power spectra, as well as the proposed scaling relation 
and relevant cosmological parameters, are the inputs to 
the cluster likelihood function. 

At this point, the analysis follows a similar path to that 
laid out by VIO. First, the matter power spectra and 
cosmology are used to calculate a mass function based 
upon the Tinker et al. (2008) prescription, which we cal- 
culate for an over-density of A = 500 flrn{z), to match 
our cluster mass definition in Section 3.1. As noted in 
Tinker et al. (2008), this function predicts the halo abun- 
dance function of input cosmology across a mass 
range of lO^/jriM© < M < IO^^/i-^Mq and a redshift 
range of < z < 2.5. Tinker et al. (2008) claim an 
overall calibration of their mass function to simulations 
of < 5%. Stanek et al. (2010) found that the inclusion 
of non-gravitational physics can shift the normalization 
of the mass function by ~ 10% along the mass direc- 
tion. However, this effect is approximately degenerate 
with an uncertainty between intra-cluster gas observables 
and mass, which we account for explicitly in our scaling 
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relation uncertainty through equations 1 and 3. 

As in VIO, the next step in the analysis is to move the 
theoretically predicted cluster abundances from their na- 
tive M500 mass space into the observable space for this 
analysis. VIO define this space by the SZ detection sig- 
nificance, f , and the optically derived rcdshift, z. This 
resulted in a two-dimensional surface of predicted cluster 
abundances in the observable space. In this analysis, we 
perform a similar transformation, this time including a 
third dimension, the X-ray parameter Yx- This results 
in a three-dimensional volume of predicted cluster abun- 
dances, now as a function of ^, Yx, and z. 

Using the scaling relations discussed in §3.1, the halo 
mass function is recast as a predicted number density in 
terms of ^, Yx and z, which we write as 

dN{i,Yx,z\v) 



dS,dYxdz 



= J dMP{tYx\M,z,p)PiM,z\p)e{^~5) 



(4) 

where p is the set of cosmological and scaling relation 
parameters, and Q is the Heaviside step function. The 
likelihood function is then given by the Poisson proba- 
bility: 



ln/:(p)-^ln 



dN{^,,Yx,,z,,\p) 



d^dYxdz 

dN{^,Yx,z,\p) 
d^dYxdz 



dS^dYxdz, 



(5) 



where the sum over the i index runs over the SPT cluster 
catalog. Note that we have neglected a global offset to 
the log-likelihood. 

We compute Equation 4 on a three-dimensional grid 
that is 200 by 200 by 30 in the C, Yx, and z dimensions, 
respectively. For each value of Yx and z we then con- 
vert to the ^ basis by using the C-^ relation defined in 
Equation 2, where we also convolve with a unit- width 
Gaussian in ^ to account for the noise in the SPT mea- 
surement. 

For each step in the MCMC, we recalculate Yx for 
each cluster given its Tx [r) and Mg (r) profiles from Sec- 
tion 2.1.3, so that its calculated Yx is consistent with 
the Yx — M500 relation and rsoo at that step. To ac- 
count for this in the cosmological likelihood, we modify 
the likelihood by adding ^ - In Yxi to the right hand side 
of Equation 5. For a detailed explanation, see Appendix 
B. For each cluster, we account for finite measurement 
errors or missing data in z and Yx by modifying the first 
term in Equation 5 by marginalizing over the relevant 
parameter, weighted by either a Gaussian likelihood de- 
termined from its uncertainty or a uniform distribution 
over the allowed range. 

From this calculation we obtain a value for the clus- 
ter likelihood corresponding to this particular set of cos- 
mological and scaling relation parameters. This value 
is then returned to CosmoMC where it may be com- 
bined with other likelihood calculations from supplemen- 
tal data sets and is used in the MCMC step accep- 
tance/rejection computation. 

4. ACDM RESULTS 

We first consider the SPTcl data constraints for a spa- 
tially flat ACDM cosmological model. For this model, we 
fit 15 parameters: the nine scaling relation parameters 
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Fig. 1. — Assuming a ACDM cosmology, the two-dimensional 
marginalized constraints on eg and ^rn- Contours show the 68% 
and 95% confidence regions for the SPTcL-t-i^o+BBN (red), CMB 
(gray), and CMB+SPTql (blue) data sets. The black lines are 
the best-fit constraint (solid) and 68% confidence region (dashed) 
for the combination of parameters that the SPTql+Hq-I-BBN 
data set best constrains: (T8(Cm/0.25)'''3° = 0.785 ± 0.037. 

and six primary cosmology parameters listed in Table 3. 
For constraints on any individual parameter, we always 
quote the mean of the likelihood distribution and the 
68% confidence interval about the mean. The confidence 
interval reflects uncertainties after marginalizing over all 
other parameters, and includes systematic uncertainties 
in the cluster scaling relations and mass calibration, as 
described in Section 3.1. In this analysis, we use the 
SPTcL and external cosmological data sets as described 
in Section 2. 

4.1. Cosmological Constraints 

The SPTcL data are not sensitive to all six ACDM 
cosmology parameters. Here and in Section 5.1, when 
considering the SPT^l cosmological constraints without 
CMB data, we always include BBN and iJo priors, as in- 
dicated. For the SPTcL+-ffo+BBN data set, we also fix 
the optical depth of reionization, r, and allow the scalar 
tilt, rig, to vary uniformly between 0.944 and 0.989, the 
95% confidence range from Keisler et al. (2011) assum- 
ing a ACDM model. However, we note that the SPTcl 
cosmological constraints vary negligibly over this range 
of Ug. As noted in Section 2.2, whenever we refer to the 
SPTcL data, we are implicitly referring to the combined 
SPT-SZ data, optical redshift, and X-ray measurements 
described in Section 2.1. 

In Figure 1, we show the constraints on the Cg and 51^ 
parameters for the individual and combined SPTcl and 
CMB data sets. In Table 3, we give the marginalized con- 
straints for the cosmological and scaling relation param- 
eters. The latter will be discussed further in Sections 4.2 
and 6. In a ACDM cosmology, the SPTcl data is most 
sensitive to cg and The number of clusters increases 
with either parameter, so the cluster abundance data ef- 
fectively constrain a product of the two. We find that 
the SPTcl +-ffo+BBN constraints are well approximated 
as a8(rJm/0.25)"-30 = 0.785 ± 0.037, which we show in 
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TABLE 2 
Parameter Table 



Type 


Symbol 


Meaning 


Gaussian Prior 


Scaling 


Asz 


C-mass normalization 


5.58 ± 1.67 


Relation 


Bsz 


C-mass slope 


1.32 ± 0.26 


Parameters 


Csz 


C-mass redshift evolution 


0.87 ±0.44 




Dsz 


Log-normal scatter in ( 


0.24 ± 0.16 




Ax 


Yx-mass normalization 


5.77 ±0.56 




Bx 


Yx-mass slope 


0.57 ±0.03 




Cx 


Yx-mass redshift evolution 


-0.40 ±0.20 




Dx 


Log-normal scatter in Yx 


0.12 ±0.08 




P 


Correlated scatter between ^ and Yx 


Uniform: (0.02, 0.98) 


Primaxy 


Qch^ 


Dark matter density 




Cosmology 




Baryon density 




Parameters 


lOOGs 


Angular scale of the sound horizon at last scattering 






Us 


Scalar tilt of power spectrum 








Scalar amplitude of power spectrum 






T 


Optical depth to reionization 




Extension 


W 


Dark energy equation of state 




Cosmology 




Fraction of dark matter in the form of neutrinos, T,mv= 94eV(/^f2c7i^) 




Parameters 




The effective number of relativistic species 






fNL 


Primordial non-Gaussianity parameter 




Derived 


o-g 


Matter fluctuations on 8 Mpc scales at z = 




Cosmology 




Total matter density 




Parameters 


h 


h = Ho/wo km s~^ Mpc~^, where Ho is the Hubble constant at z = 





TABLE 3 
ACDM Constraints 



Parameter 


Prior 


SPTcL+i?o+BBN 


CMB 


CMB+SPTcL 


Asz 


5.58 ± 1.67 


5.31 ± 0.98 




4.91 ± 0.71 


Bsz 


1.32 ±0.26 


1.39 ±0.15 




1.40 ±0.15 


Csz 


0.87 ± 0.44 


0.90 ±0.34 




0.83 ±0.30 


Dsz 


0.24 ±0.16 


0.21 ±0.10 




0.21 ±0.09 


Ax 


5.77 ±0.56 


5.69 ±0.51 




5.82 ± 0.48 


Bx 


0.57 ±0.03 


0.564 ±0.029 




0.563 ±0.029 


Cx 


-0.40 ±0.20 


-0.37 ±0.16 




-0.35 ±0.16 


Dx 


0.12 ±0.08 


0.14 ±0.07 




0.14 ±0.07 


P 


(0.02,0.98) 


0.52 ±0.27 




0.52 ±0.27 






0.133 ± 0.045 


0.111 ±0.0048 


0.109 ±0.0032 






0.0221 ±0.0020 


0.0222 ± 0.0004 


0.0223 ± 0.0004 


looe^ 




1.065 ± 0.041 


1.041 ±0.0016 


1.041 ±0.0016 


Us 


(0.944, 0.989) 


0.966 ± 0.013 


0.965 ± 0.011 


0.967 ±0.010 


lO^Als 




2.16 ± 1.30 


2.44 ±0.10 


2.40 ± 0.08 


r 




(0.090) 


0.086 ±0.014 


0.087 ±0.014 






0.766 ± 0.062 


0.808 ± 0.024 


0.795 ±0.016 






0.285 ± 0.083 


0.268 ±0.025 


0.255 ±0.016 


h 




0.739 ± 0.024 


0.707 ±0.022 


0.717 ±0.016 



Note. — The marginalized constraints on the scaling relation and primary cosmology parameters from Table 2, where we 
report the mean of the likelihood distribution and the 68% confidence interval about the mean. The priors are Gaussian, except 
for p and Us, which are uniform over the range given. The Us prior is only used for the SPTcL±^fo±BBN data set. 
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[10^^M3,„ keV] 

Fig. 2. — A plot of the SZ-significance, t;, versus the X-ray 
observable Yx for the 14 SPT clusters with X-ray measurements. 
From the CMB-t-SPTpL data set fit to a ACDM cosmology, 
we use the best-fit f — M500 and Yx — M500 scaling relations 
to calculate the expected form and redshift evolution of the 
^-Yxrelation (solid-line), where Epiv = E{z = 0.6). With the 
best-fit cosmology parameters, we also predict the effective 68 and 
95% confidence intervals for the expected distribution of clusters 
in the ^-Yx plane (red contours). The measured and predicted 
cluster distribution show qualitatively good agreement. 

Figure 1 by the solid and dashed Hnes. Combining the 
SPTcL and CMB data, we constrain erg = 0.795 ± 0.016 
and flm = 0.255 ± 0.016, a factor of 1.5 improvement on 
each over the constraints from the CMB alone. 

The SPTcL constraints are consistent with results us- 
ing optical and X-ray selected cluster samples. Re- 
cently, Rozo et al. (2010) compared the cluster con- 
straints from several different methods, and found gen- 
erally good agreement and comparable constraints. It 
is typical for cluster based constraints to be quoted in 
terms of the product of erg and Qm to an exponent which 
varies depending on the mass scale of the cluster sample. 
One example for comparison is Vikhlinin et al. (2009b), 
who constrained cr8(17„/0.25)° ''^ = 0.813 ± 0.027. For 
typical ACDM model constraints of ^ 0.25 — 0.30, 
this agrees well with our result. 

4.2. Scaling Relation Constraints 

In Figure 2, we show the relationship between ^ and Yx 
for the 14 clusters with X-ray observations, over-plotted 
with the expected distribution of clusters and the best-fit 
relation determined using the CMB+SPTcl data from 
Section 4.1. The combination of the steep mass function 
and SZ selection yields a distribution of clusters visibly 
offset from the best-fit scaling relation, an effect often re- 
ferred to as Eddington bias. We note that our cosmolog- 
ical analysis method described in Section 3.2, explicitly 
accounts for the SZ selection and therefore Eddington 
bias. We also predict the expected distribution of clus- 
ters in the S,-Yx plane assuming the best-fit cosmology 
and scaling relation parameters, and applying a compa- 
rable selection as was used for the SPT X-ray follow-up 
(z > 0.3 and ^> 5.45). The predicted 14.2 clusters is 



consistent with the 14 detected. In Figure 2, we over- 
plot the effective 68 and 95% confidence region in the 
^-Yx plane where we would expect to find these clus- 
ters. Qualitatively we find good agreement between the 
observed and predicted cluster distribution. 

In Table 3, we give the constraints on the Yx ~ A/500 
and C — -^'^500 scaling relations using the SPTcl and 
CMB+SPTcL data sets. Because the Yx - A'hoo rela- 
tion has significantly tighter priors than the ( — M^qq 
relation, we will not give the Yx — M500 constraints for 
the modified cosmologies presented in Section 5. Simi- 
larly, for the parameter p, the correlated scatter between 
C and Yx, we have virtually no constraining power. In 
all cases, p moves nearly uniformly across the entire al- 
lowed range, and has a negligible effect on the cosmolog- 
ical constraints. In Appendix C, we give posterior mass 
estimates for each cluster using a similar method as de- 
scribed in VIO and briefly reviewed in the appendix. 

If there were a significant discrepancy between the 
simulation-based prior on the ^ — M500 relation and 
the observational prior on the Yx — M^oo relation, we 
would observe it as an offset between the central value 
of the C — -M500 prior and its best-fit value. From the 
CMB+SPTcL constraints, the largest offset is for Agz, 
with a best-fit value of 4.91 ± 0.71 compared to the sim- 
ulation prior of 5.58 ± 1.67. An offset in this direc- 
tion would be consistent with the SZ simulation prior 
under-estimating the mass of a cluster by a factor of 
- ((4.91 ± 0.71)/5. 58)1/1'' ^ Q 9;^ ^ 0.09. This result is 
consistent with preliminary estimates from All, who es- 
timated this factor to be 0.78 ± 0.06. We note that All 
did not marginalize over uncertainties in either the X- 
ray scaling relation or cosmological parameters, both of 
which affect this result. The derived offset is also a func- 
tion of the assumed cosmology. For example, if we as- 
sume a ACDM cosmology with a non-zero neutrino mass, 
as in Section 5.2, we find a value closer to the simulation 
prior, Asz = 5.39 ± 0.79, using the CMB-^Ho+SPTcl 
data set. Therefore, we find no significant inconsistency 
between the simulation-based prior on the C — M500 re- 
lation and the observational prior on the Yx — M^qq re- 
lation. 

5. EXTENSIONS TO ACDM 

In this section, we consider extensions to a spatially 
flat ACDM cosmology. For each extension, we also fit the 
nine scaling relation parameters and six primary cosmol- 
ogy parameters listed in Table 2. We consider four exten- 
sion cosmologies where we include the following as free 
parameters: the dark energy equation of state (w), the 
sum of the neutrino masses (Em^), the sum of neutrino 
masses and the effective number of relativistic species 
(A'cff ), and a primordial non-Gaussianity {/nl)- For con- 
straints on any individual parameter, we always quote 
the mean of the likelihood distribution and the 68% con- 
fidence interval about the mean. The confidence inter- 
val will include uncertainties after marginalizing over all 
other parameters, which includes systematic uncertain- 
ties in the cluster scaling relations and mass calibration, 
as described in Section 3.1. In this analysis, we use the 
SPTcL and external cosmological data sets as described 
in Section 2. 

5.1. wCDM 
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TABLE 4 
wCBM Constraints 







CMB 


Hq+BBN 


CMB+BAO+SNe 


CMB+BAO+SNe 








+SPTCL 




+SPTCL 


Scaling 


Asz 




5.12 ± 1.36 




4.75 ± 0.79 


Parameters 


Bsz 




1.40 ±0.15 




1.41 ± 0.15 




Csz 




0.92 ±0.36 




0.85 ± 0.29 




Dsz 




0.22 ±0.10 




0.21 ± 0.10 


Cosmology 




0.864 ±0.120 


0.773 ±0.088 


0.823 ±0.040 


0.793 ±0.028 


Parameters 




0.244 ±0.089 


0.293 ±0.113 


0.279 ±0.016 


0.273 ±0.015 




h 


0.775 ±0.128 


0.740 ± 0.025 


0.698 ± 0.018 


0.697 ±0.018 




w 


-1.19 ±0.37 


-1.09 ±0.36 


-1.014 ±0.078 


-0.973 ± 0.063 



Note. — The marginalized constraints on a subset of the scaling relation and cosmology parameters from Table 2. Scaling relation and 
primary cosmology parameters not given are still varied in the MCMC and marginalized over for these constraints. We report the mean of 
the likelihood distribution and the 68% confidence interval about the mean. 



TABLE 5 
ACDM + 'Em„ + iVefi Constraints 







CMB+Ho+BAO 


CMB+Ho+BAO 


CMB+Ho 


CMB+Ho+BAO 








+SPTCL 


+SPTcL 


+SPTcL 


Scaling 


Asz 




5.26 + 0.79 


5.39 ± 0.79 


5.01 ± 0.85 


Parameters 


Bsz 




1.39 ±0.14 


1.39 ±0.14 


1.41 ± 0.15 




Csz 




0.89 ±0.30 


0.89 ±0.30 


0.89 ±0.30 




Dsz 




0.20 ±0.09 


0.20 ±0.10 


0.21 ± 0.10 


Cosmology 




0.761 ± 0.043 


0.770 ±0.026 


0.771 ± 0.023 


0.777 + 0.031 


Parameters 




0.275 ±0.016 


0.272 ±0.015 


0.260 ±0.018 


0.284 + 0.018 




h 


0.698 ±0.014 


0.701 ± 0.013 


0.712 ± 0.017 


0.727 + 0.020 






0.19 ±0.14 


0.15 ±0.10 


0.12 ±0.09 


0.34 ±0.17 




T,m^(eV), 95% CL 


< 0.45 


< 0.33 


< 0.28 


< 0.63 






(3.046) 


(3.046) 


(3.046) 


3.91 ± 0.42 



Note. — The marginalized constraints on a subset of the scaling relation and cosmology parameters from Table 2. Scaling relation and 
primary cosmology parameters not given are still varied in the MCMC and marginalized over for these constraints. We report the mean of 
the likelihood distribution and the 68% confidence interval about the mean. 
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Fig. 3. — Assuming a u;CDM cosmology, the two-dimensional 
marginalized constraints on w and erg. Contours show the 68% 
and 95% confidence regions for the SPTgL+Z/o+BBN (red) and 
CMB (gray) data sets. 

The first extension we consider is a wCDM cosmology, 
a model in which the equation of state of dark energy is 
a constant w. The cluster abundance and the shape of 
the mass function depend on w through its effect on the 
growth of structure, or cquivalcntly the redshift evolution 
of (Tg. The CMB data measure structure at z ^ 1100, 
and therefore require significant extrapolation to predict 
the cluster abundance in the redshift range of the SPT 
sample (0.3 < z < 1.1). Therefore, consistency between 
the implied w from both data sets is already an important 
systematic test of dark energy. 

In Figure 3, we show the constraints on w and erg using 
the CMB and SPTcL+i?o+BBN data sets. The likeli- 
hood contours have significant overlap, implying the data 
are in good agreement. Relative to the CMB, the SPTcl 
data tend to disfavor cosmologies with large erg and more 
negative w. In Table 4, we give marginalized constraints 
for several cosmological and scaling relation parameters. 
The SPTcL data constrain w = —1.09 ± 0.36, and have 
similar constraining power to the CMB data, for which 
the constraints have a significant degeneracy between 
w and (Tg. The SPT^l data simultaneously constrain 
(Tg = 0.773 ± 0.088. This constraint has a factor of ^lA 
lower uncertainty than that from the CMB data. 

5.1.1. wCDM with BAO and SNe data sets 

In this section, we consider the improvement in wCDM 
cosmological constraints when adding the SPTcl data to 
the CMB, BAO, and SNe data sets. In Figure 4, we show 
the constraints of the combined CMB+BAO-|-SNe data 
set, before and after including the SPTcl data. The 
SPTcL data most significantly improve the constraints 
on (Tg and w; reducing the allowed two-dimensional like- 
lihood area by a factor of ~1.8. In Table 4, we give 
the marginalized constraints for several parameters be- 
fore and after the inclusion of the SPTcl data. The 
combined constraints are w = —0.973 ± 0.063 and erg = 
0.793 ± 0.028, a factor of 1.25 and 1.4 improvement, re- 
spectively, over the constraints without clusters. The 



combined data set also constrains Q.m = 0.273 ± 0.015 
and h — 0.697 ± 0.018. These constraints are consis- 
tent with previous cluster-based results (Vikhlinin et al. 
2009b; Mantz et al. 2010c; Rozo et al. 2010), which used 
X-ray and optically selected samples of typically lower 
redshift clusters. The sensitivity of the SPTcl clus- 
ter data to the amplitude of structure, erg, is primarily 
what gives it the ability to break degeneracies with the 
distance-relation based constraints from the BAO and 
SNe data sets. We note the slight tension with the Hq 
constraints from Riess et al. (2011) of h = 0.738 ± 0.024. 
While this tension is not significant, it helps to intuitively 
explain some constraints on neutrino mass in Section 5.2. 

5.1.2. — A/500 Constraints 

Given the work of VIO and other cluster results (e.g., 
Vikhlinin et al. 2009b; Mantz et al. 2010c; Rozo et al. 
2010), we expect the cluster mass-calibration to be the 
dominant systematic uncertainty limiting our results. In 
Figure 5, we show the constraints on Asz and erg. The 
SPTcl +-f^o+BBN data set has a significant degeneracy 
between its constraints on A^z and erg. From this data, 
we constrain the fractional uncertainty, 5 Asz I -^szi to be 
27%, which is effectively constrained only by the 14 clus- 
ters that have both X-ray and SZ measurements. This 
constraint is not significantly better than the uncertainty 
in the simulation based prior of 30%. With enough X- 
ray observations, we expect the — M500 calibration to 
be limited by the uncertainty of the Yx — M^qq relation, 
because the latter is currently better observationally con- 
strained. In this limit, we would expect a fractional un- 
certainty on Asz of Bsz X {SAx I Ax) ~ 14%. The above 
would suggest that for a wCDM cosmology we would 
need X-ray observations of ^ 50 clusters, i.e., ~ 14 clus- 
ters x(27%/14%)^, to calibrate Asz in terms of mass so 
that its not the dominant source of uncertainty. When 
adding the BAO and SNe data sets, we improve the con- 
straints on Asz to an accuracy of ^ 16%. However, these 
data sets are not sensitive to either Asz or erg, and can- 
not completely break their degeneracy. In Section 6, we 
will discuss the systematic uncertainties from this degen- 
eracy on our cosmological constraints in more detail. 

5.2. KCDM with Massive Neutrinos 

We next consider a ACDM cosmology with non-zero 
neutrino masses. Cosmological measurements are pri- 
marily sensitive to the neutrino masses through their ef- 
fect on structure formation. A massive neutrino addi- 
tionally affects the CMB power spectrum if it was non- 
relativistic at the redshift of recombination. For exam- 
ple, if the heaviest neutrino had a mass < 0.6 eV, it would 
be relativistic at recombination and therefore would not 
significantly affect the structure in the CMB (Komatsu 
et al. 2009). However, as the universe expanded and 
cooled neutrinos would transition to non-relativistic, and 
would contribute to but not to structure formation 
below their free streaming scale, implying a lower erg at 
z — and fewer clusters. This implies that measure- 
ments of the CMB power spectrum alone cannot con- 
strain the neutrino mass to significantly less than 0.6 eV 
per species (i.e., Sm^< 1.8 eV), and the constraints will 
be significantly degenerate with erg. Local measurements 
of structure break this degeneracy, and significantly im- 
prove the neutrino mass constraints. 
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Fig. 4. — Assuming a uiCDM cosmology, the constraints on Qm, eg, and w. The plots along the diagonal are the one-dimensional 
marginalized likelihood. The off-diagonal plots are the two-dimensional marginalized constraints showing the 68% and 95% confidence 
regions. We show the constraints for the CMB-|-BAO-f SNe (gray, dashed), and CMB+BAO-f-SNe-f SPTql (green, solid) data sets. The 
SPTql data improves the constraints on erg and w, by factors of 1.4 and 1.25, respectively. 



We also note the significant degeneracy between the 
CMB power spectrum constraints on Hq and Em^ . Mas- 
sive neutrinos affect the amplitude of the early integrated 
Sachs- Wolfe effect causing a shift of the first peak of 
the CMB power spectrum towards larger angular scales 
that can be absorbed by a lower value of Hq (Ichikawa 
et al. 2005). Komatsu et al. (2011) used a combination 
of WMAP7-t-i7o+BAO data to set a hmit of i:m^< 0.58 
at a 95% confidence limit (CL). Following Komatsu et al. 
(2011), we consider the same combination of data sets to 
add to the CMB power spectrum measurements, which 
were chosen because of their insensitivity to systematic 
errors and their ability to maximally constrain Em^ by 
breaking the degeneracy with Hq. We define 'Sm^— 94 
eV{f^il.ch^), where fi, is the fraction of dark matter in 
the form of massive neutrinos. 

In Figure 6, we show the constraints on cig and Etoi^, 
using the CMB+_ffo+BAO data set, before and after 
including the SPTcl data. In Table 5, we give the 
marginalized constraints on each parameter. Using the 



CMB-fHo+BAO+SPTcL data set, we constrain T,mi,< 
0.33 eV at a 95% CL, a factor of 1.4 improvement over 
the constraints without the SPTcl data. This improve- 
ment is primarily due to the tighter constraints on as for 
which the uncertainty decreased by a factor 1.8. The con- 
straint is lower by excluding the BAO data; using only 
the CMB+i/o+SPTcL data set we constrain Em^< 0.28 
eV at a 95% CL. These improved constraints can be un- 
derstood from the Hq measurements, as also noted in 
Section 5.1. The results of Riess et al. (2011) favor a 
marginally higher Hq value than the CMI3+BA0 data. 
Because of the degeneracy between Em^ and Hq in the 
CMB constraints, a higher value of Hq tends to favor 
lower values of 'Sm^. The constraints presented here are 
comparable to other recent results using optically and 
X-ray selected cluster samples with similar cosmological 
data sets (Reid et al. 2010; Mantz et al. 2010b). 

5.2.1. Number of Relativistic Species 
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Fig. 5. — Assuming a u;CDM cosmology, the two-dimensional 
marginalized constraints on Agz a^nd ""S- Contours show the 
68% and 95% confidence regions for the SPTcl+^^O+BBN (red) 
and CMB+BAO+SNe+SPTcL (green) data sets. The horizontal 
black dashed line is the center of the theory prior on Agz- 
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Fig. 6. — Assuming a ACDM cosmology with massive neutrinos, 
the two-dimensional marginalized constraints on Em^ and 
erg. Contours show the 68% and 95% confidence regions for 
the CMB-I-Ho+BAO (gray, dashed), CMB-|-Ho-|-BAO-(-SPTcl 
(orange, solid), and CMB-|-Ho-|-SPTcl (blue, dot-dashed) data 
sets. The SPT(^l data improves the constraints on ag and Sm^, 
by factors of 1.8 and 1.4, respectively. 

Recent measurements have shown a ~2a preference for 
increased damping in the tail of the CMB power spec- 
trum (Dunkley et al. 2011; Keisler et al. 2011). This 
damping could be caused by several different physical 
mechanisms, such as a high primordial helium abun- 
dance, a running of the scalar spectral index, or addi- 
tional relativistic species. This last explanation is partic- 
ularly timely because of recent measurements from atmo- 
spheric (Aguilar-Arevalo et al. 2010) and nuclear reactor 
(Mention et al. 2011) neutrino oscillation experiments 



that find some evidence for a sterile neutrino species. It 
has been pointed out that these measurements are most 
consistent with two sterile neutrinos and Stoi/> 1.7 eV 
(Kopp et al. 2011). Therefore, we consider the joint cos- 
mological constraints on iVoff and Y,m^ to compare with 
these terrestrial results. 

With only three neutrino species, we would expect 
Neff= 3.046, a value slightly larger than three because 
of energy injection from electron-positron annihilation 
at the end of neutrino freeze-out (Dicus et al. 1982; 
Lopez et al. 1999; Mangano et al. 2005). As TVoff in- 
creases, the contribution to the gravitational potential 
of the additional neutrino perturbations boosts the early 
growth of dark matter perturbations (Bashinsky & Sel- 
jak 2004), which also increases as (Hou et al. 2011). As 
explained in Section 5.2, adding neutrino mass at the 
levels considered here only affects the low-redshift uni- 
verse, suppressing structure formation, and lowering erg 
at z = 0. Therefore, increasing N^ff will also allow an in- 
creasing 'Em^. Keisler et al. (2011) used a combination 
of CMB+i?o+BAO data to constrain Em^< 0.69 eV at 
a 95% CL, erg = 0.803 ± 0.056, and TVeff = 3.98 ± 0.43. 

In Figure 7, we show the constraints on N^g^ Smi,, 
and CTg, using the CMB-fi/o+BAO data set, before and 
after including the SPTcl data. In Table, 5 we give 
the marginalized constraints. When varying A^eff we as- 
sume consistency with BBN for our constraints. Us- 
ing the CMB-l-i/o+BAO+SPTcL data set, we constrain 
I]m^< 0.63 eV at a 95% CL, erg = 0.777 ± 0.031, and 
A^,ff= 3.91 ± 0.42. Relative to Keisler et al. (2011), the 
addition of the SPTcl data improves the constraints on 
cTg by a factor of 1.8, and reduces the upper limit on T,mi, 
by a factor of 1.1. However, the addition of the SPTql 
data does noticeably sharpen the peak in the marginal- 
ized one-dimensional likelihood for 'Sm^, such that the 
maximum likelihood constraint peaks away from zero, 
I]mj,= 0.34 ±0.17 eV. 

As noted in Keisler et al. (2011), models of the CMB 
power spectrum that include increased damping are fa- 
vored at the 1.6-1.9(7 level. However, even if one accepts 
the need for an extra parameter to explain the damping, 
its physical origin is unclear. Regardless, considering the 
A'cff model extension is instructive to help understand 
the model dependency of the neutrino mass constraints. 
Keisler et al. (2011) considered three models to explain 
the excess damping and found that the N^s model had 
the most significant effect on dg. The inclusion of N^s 
also weakens the constraints on Smj^, because of the de- 
generacies between A^eff, cg, and Sm^. In the combined 
cosmological data set, the SPTcl data mainly constrains 
(Tg, which helps to break this degeneracy and indirectly 
improve the neutrino mass constraints. Therefore, the 
Emi^ constraint from the Nqs model can be considered a 
conservative upper limit on Sm^ regardless of the phys- 
ical mechanism for the increased damping. 

5.3. ACDM with Jnl 

Finally, we consider a ACDM cosmology with primor- 
dial non-Gaussianity. Standard inflationary cosmology 
predicts that density fluctuations in the universe were 
seeded by random Gaussian fluctuations. However, in- 
flationary models can be constructed that predict sig- 
nificant levels of non-Gaussianity (e.g., Bartolo et al. 
2004). The leading order non-Gaussian term is typi- 
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Fig. 7. — Assuming a ACDM cosmology with massive neutrinos and a free number of relativistic species, N^g, the constraints on trg, Sm^, 
and N^ff. The plots along the diagonal are the one-dimensional marginalized likelihood. The off-diagonal plots are the two-dimensional 
marginalized constraints showing the 68% and 95% confidence regions. We show the constraints for the CMB-f Hq-I-BAO (gray, dashed), and 
CMB-(-Ho-|-BAO-f SPTql (orange, solid) data sets. The black vertical line shows A^eff = 3.046, the expected value for three neutrino species. 
Using the CMB+Ho-f BAO-(-SPTcl data set, the 68% confidence marginalized constraints are erg = 0.777 it 0.031, A^cff = 3.91 it 0.42, and 
Sm^= 0.34 it 0.17 eV, with a 95% CL of T,m^< 0.63 eV. 



cally described by the parameter /atl- Using a mea- 
surement of the CMB power spectrum from WMAP7 
data, Komatsu et al. (2011) measured a 95% CL of 
— 10 < Jnl < 74. Primordial non-Gaussianity can also 
manifest itself through the abundance of massive galaxy 
clusters. While the constraints from an SPT-like SZ sur- 
vey are not expected to be competitive with current CMB 
constraints (Dalai et al. 2008), they constitute a con- 
straint independent from the CMB results which is sen- 
sitive to very different physical scales. In principle, even 
a single massive high-redshift cluster can falsify standard 
ACDM cosmology (Mortonson et al. 2011), though cur- 
rently the most massive cluster known at z > 1 is not in 
significant tension (Foley et al. 2011). In our analysis, 
we incorporate /jv^ as a modification of the cluster mass 
function following the prescription of Dalai et al. (2008), 
and only consider the effect of /nl on the SPTcl data 
set. 

In Figure 8, we show the constraints on erg and /nl 



using the CMB-|-SPTcl data set. The marginalized con- 
straints are fNL= -192 ± 310, as = 0.803 ± 0.022, and 
Asz = 5.27 ± 0.89. The best-fit value of Jnl is slightly 
negative, generally implying fewer massive high-redshift 
clusters. In Section 4, we found that for a ACDM cos- 
mology using the CMB+SPT^l data, our constraint for 
Asz was 0.9 a lower than the simulation prior. This cor- 
responds to a higher mass for a cluster of a given ^, which 
can also be thought of as the SPT survey having fewer 
clusters than expected given the simulation prior. When 
/jVL is added as a parameter, Asz moves back towards 
its simulation prior, and the deficit of clusters can be 
maintained by a more negative Jnl- For the model to 
match the number of clusters measured by SPT, a more 
negative /nl can be balanced with either a larger erg or 
Asz, creating a degeneracy in this direction. Regard- 
less, any deficit of clusters is not significant relative to 
the uncertainty on either erg or Asz, even in a ACDM 
cosmology for which they are best constrained. 
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Fig. 8. — Assuming a ACDM cosmology with a primordial 
non-Gaussianity characterized by the parameter /]vli the two- 
dimensional marginalized constraints on erg and /jvl using the 
CMB+SPTcL data set (blue). Contours show the 68% and 95% 
confidence regions. We only consider the affect of /jvl on the 
SPT(^L data set. We measure /jvl= —192 it 310, consistent with 
zero non-Gaussianity. 

Our results are consistent with Williamson et al. (2011) 
which used the 26 most massive clusters in the full 2500 
deg2 SPT-SZ survey to constrain fNL= 20 ± 450. Our 
work differs from Williamson et al. (2011) in that we use 
a much smaller area of the SPT-SZ survey, we select clus- 
ters down to a lower ^ threshold, and we use an improved 
mass calibration. 

6. SOURCES OF UNCERTAINTY 

Previous SPT cluster survey results, namely VIO and 
Williamson et al. (2011), found cosmological constraints 
that were limited most significantly by the cluster mass 
calibration, or equivalently the fractional uncertainty in 
Asz- In this work, we have reduced this uncertainty 
by incorporating the external mass calibration from the 
Yx ~ M500 relation using X-ray observations of the SPT 
clusters. We can directly estimate the impact of the 
uncertainties in the X-ray and SZ scaling relations by 
importance sampling the MCMC chains, where we post- 
process the chains by imposing a narrow prior on each 
scaling relation parameter centered around the best-fit 
value. The resulting increase in precision on the cosmo- 
logical parameters allows a measure of the impact from 
the uncertainty in the scaling relations. In this way, we 
effectively "fix" the X-ray and scaling relation param- 
eters, a process which we will implicitly be referring to 
throughout this section. For a wCDM cosmology, we also 
consider the impact of the SNe systematic uncertainty on 
the cosmological results presented here. 

With enough SZ and X-ray observations, we expect 
the C — -^^500 calibration to be limited by the calibra- 
tion of the Yx — M500 relation because the latter has 
tighter external priors. In practice, there will be an ad- 
ditional uncertainty in the ^ — M500 calibration from the 
limited number of SZ and X-ray observations for cross- 
calibration, and this uncertainty will also degrade the 
cosmological constraints. We wish to separate this ef- 



fect, which we will refer to as the SZ-Yx scaling uncer- 
tainty, from the additional systematic uncertainty from 
the Yx — M500 calibration, which we will refer to as the 
X-ray scaling uncertainty, and the statistical uncertainty 
from the cluster sample size. By fixing the X-ray and SZ 
scaling relation parameters, as described above, we can 
measure the impact of the SZ-Yx scaling uncertainty, X- 
ray scaling uncertainty, and statistical uncertainty on our 
cosmological constraints. 

6.1. ACDM Cosmology: Scaling Relation Uncertainty 

We first consider the ACDM constraints using the 
SPTcL+-f^o+BBN data set, the results of which were 
described in Section 4. This data best constrained the 
combination of cr8(fi.,„/0.25)°-30 = 0.785 ± 0.037. The 
sources of uncertainty for this result are summarized in 
Table 6, and are discussed below. 

For the X-ray scaling relation parameters, only the un- 
certainty in Ax and C'x, the normalization and redshift 
evolution parameters, contribute significantly to the un- 
certainty on (Tg(17m/0.25)°'^°. Fixing each parameter 
separately implies that they contribute an uncertainty 
on cr8(f^m/0.25)°-30 of ±0.022 and ±0.015, respectively. 
It is not surprising that the normalization of the mass 
calibration significantly affects the constraints, and the 
redshift evolution can be understood for similar reasons. 
For a cluster at the median redshift of the SPT sam- 
ple, z — 0.74, the prior on the Cx value effectively con- 
tributes an additional 7% to the cluster mass calibra- 
tion. This can be compared to the 10% mass calibration 
uncertainty from the prior on Ax- Fixing all X-ray pa- 
rameters simultaneously, implies that they contribute an 
uncertainty on cr8(f^m/0.25)°-3° of ±0.028. 

For the SZ scaling relation parameters, only the uncer- 
tainty in Asz contributes significantly to the uncertainty 
on (T8(f^m/0.25)°'^''. Fixing all the SZ scaling parame- 
ters, we measure an uncertainty on cr8(^m/0.25)° '^° of 
±0.023 from statistical uncertainty, ±0.010 from the SZ- 
Yx scaling uncertainty, and ±0.028 due to X-ray scal- 
ing uncertainty (as discussed above). The relatively low 
contribution from the SZ-Yx scaling uncertainty is not 
surprising considering the constraints on the fractional 
uncertainty of Agzj which was near the systematic limit 
of 14% imposed by the Yx — M^qq calibration. 

Therefore, the ACDM constraints are nearly at the sys- 
tematic limit from the calibration of the Yx — M^qq re- 
lation. For our constraint of cr8(^^m/0.25)°-^° = 0.785 ± 
0.037, the X-ray scaling and statistical uncertainty con- 
tribute almost equal amounts of ±0.028 and ±0.023, re- 
spectively. By only increasing the cluster sample size 
we could reduce the uncertainty by up to ^ 1.3 ('-^ 
0.037/0.028). Further improvements would require a 
more accurate cluster mass calibration. 

6.2. wCDM Cosmology: Scaling Relation and SNe 
Uncertainty 

We next consider the sources of uncertainty for 
the wCDM cosmology discussed in Section 5.1. We 
will concentrate on using the SPTcl±^o+BBN and 
CMB+BAO+SNe+SPTcL data sets, which produce con- 
straints oiw = -1.09 ± 0.36 and w = -0.973 ± 0.063, 
respectively. The sources of uncertainty for this result 
are summarized in Table 6, and are discussed below. 
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TABLE 6 
Error Budget 





ACDM 


wCDM 




iL'CDM 




SPTCL+ 


Ho+BBN 


SPTcL±^^o±BBN 


CMB+BAO+ 


SNc+SPTcL 






/0.25)0-30 


w 




w 


OS 


Baseline, Section 4,5.1 


0.785 ± 0.037 


-1.09 ±0.36 0.773 ±0.088 


-0.973 ± 0.063 


0.793 ± 0.028 


SNe Systematic 










±0.026 


±0.005 


SZ-Yx Scaling 




±0.010 


±0.19 


±0.066 


±0.013 


±0.013 


X-ray Scaling Systematic 




±0.028 


±0.15 


±0.036 


±0.019 


±0.014 


Statistical 




±0.023 


±0.27 


±0.046 


±0.033 


±0.013 



Note. — We give the mean and 68% confidence intervals for a subset of the cosmological parameters reported in Sections 4 and 5.1. 
The last four rows give the Icr error in each cosmological parameter due to the stated uncertainty. 



We first consider the SPTcL+-ffo+BBN data set. For 
the X-ray scaling relation parameters, we again find that 
the uncertainty in Ax and Cx contribute the largest un- 
certainty on w. Fixing each parameter independently 
implies they contribute an uncertainty of = ±0.10 
and ±0.11, respectively, and a total uncertainty of ±0.15. 
For the SZ scaling relation parameters, only the uncer- 
tainty in Asz contributes significantly to the uncertainty 
on w. Fixing all the SZ scaling parameters, we measure 
an uncertainty on w of: ±0.27 from statistical uncer- 
tainty, ±0.19 from the SZ-Yx scaling uncertainty, and 
±0.15 due to X-ray scaling uncertainty. Therefore, un- 
like the ACDM case, we find that our constraints on w 
would be significantly improved by adding more clusters 
and additional Yx measurements. A similar conclusion 
is reached repeating the above analysis for erg. In prin- 
ciple, adding more clusters and Yx measurements would 
reduce the uncertainty on w and as to values limited by 
the X-ray scaling uncertainty. In this limit we should 
measure w and us with an unceratinty of Sw = ±0.15 
and Sag = ±0.036, or ~ 2.5 times better than our cur- 
rent constraints. 

When considering the CMB+BAO+SNe+SPTcL data 
set. we reach qualitatively similar conclusions, however 
the total uncertainty is significantly lower because of the 
parameter degeneracies that are broken from the addi- 
tional data sets. We first re-run the MCMC chains with- 
out SNe systematic uncertainty. Fixing all the SZ scal- 
ing parameters, we measure a statistical uncertainty of 
6w = ±0.033, a factor of two improvement relative to the 
constraints including all systematic uncertainties. Com- 
paring this uncertainty to that with no fixed parame- 
ters, we estimate an uncertainty on w of ±0.019 from X- 
ray scaling tmcertainty and ±0.013 from SZ-Yx scaling 
uncertainty. The addition of the SPTcl data also sig- 
nificantly reduces the systematic uncertainty from SNe. 
Running a CMB+BAO+SNe MCMC chain with and 
without SNe systcmatics, we measure w = — 1.014±0.078 
and w = —1.017 ± 0.050, respectively. This implies 
that SNe systematics are contributing an uncertainty of 
Sw = ±0.060. After adding the SPTcl data the uncer- 
tainty from SNe systematics is reduced to 6w = ±0.026, 
a factor of 2.3 improvement. 

6.3. Point Source Contamination 

In VIO, it was argued that point source contamination 
contributed a negligible level of uncertainty relative to 
the statistical precision of the cluster sample. Since we 
are using the same cluster sample, we expect the same 
conclusion to hold, though we briefly summarize their 



arguments here. From Poisson distributed sources, the 
probability of a chance superposition of a bright point 
source (> 6 mJy) with a cluster is negligible, given the 
sky density of sources at 150 GHz (~ 1 deg~^, Vieira 
et al. (2010)). Furthermore, our cosmological analy- 
sis in Section 3 explicitly accounts for a Poisson dis- 
tributed background of sources, and the X-ray measure- 
ments gives an additional systematic check on an offset 
in the SZ measurements. Correlated emission from clus- 
ter members could potentially fill in cluster decrements. 
However, correlated radio emission has previously been 
calculated to be negligible at 150 GHz for clusters of the 
typical SPT mass scale and redshift range (Lin et al. 
2009; Schgal et al. 2010). In VIO, it was also argued 
that the level of correlated dusty emission is negligible, 
from the known quenching of star formation in massive 
clusters (Hashimoto et al. 1998), and the sub- millimeter 
luminosity function (Pascalc et al. 2009). In addition, 
more recent Spitzcr infrared observations of a sample of 
X-ray selected groups and low-mass clusters, found that 
correlated dusty emission is insignificant compared to the 
SZ signal (George et al. 2011). These arguments apply 
to the clusters in the SPT sample, which span a similar 
redshift range and are of higher mass than the objects 
considered in that work. 

7. DISCUSSION 

7.1. Improvement Relative to VIO 

In this work, the cluster sample is the same as used in 
VIO for their cosmological analysis. However, we have 
improved the cosmological constraints relative to VIO 
by including X-ray measurements in order to reduce the 
cluster mass calibration uncertainty. It is not straightfor- 
ward to quantify the improvement for two main reasons. 
First, VIO used pre-existing WMAP7 MCMC chains 
from Komatsu et al. (2011) that they importance sam- 
pled by re-weighting the chains by the likelihood of the 
SPT cluster catalog given each set of parameters. In this 
work, we generated new MCMC chains while simulta- 
neously fitting both data sets. Second, each result uses 
somewhat different external data sets, in particular for 
the CMB power spectrum measurements, where the re- 
sults in this work also include CMB measurements from 
Keisler et al. (2011). 

Without explicitly correcting for these differences, we 
can approximate the improvement from including the X- 
ray measurements by considering the relative improve- 
ments of adding the SPTcl data to the CMB data used 
in either result. For a ACDM cosmology, the WMAP7 
data constrains as = 0.801 ± 0.030. In VIO, the ad- 
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dition of the SPT data, this constraint improved to 
as = 0.791 ± 0.027, a factor of 1.1 improvement. In 
this work, for a ACDM cosmology using the SPTcl and 
CMB data, we constrained as = 0.795 ± 0.016, a factor 
of 1.5 improvement over the constraints from the CMB 
alone. Therefore, the addition of the X-ray measure- 
ments improved tlie ACDM constraints on erg by a factor 
of ^1.4. A comparison of the wCDM cosmological con- 
straints is more complicated because of the somewhat 
different handling of the external data sets. Regardless, 
the significant improvement in the constraints from the 
X-ray measurements is clear. 

7.2. Prospects for Further Improvement 

The results in this paper were derived using 18 clus- 
ters from 178 deg^ of the 2500 deg^ SPT-SZ survey. 
The full survey will significantly increasing both the 
area and overall depth of the SZ maps. Reichardt et 
al. (2012, in prep.) will present a catalog of ^200 
clusters from the first 800 deg^ of the SPT-SZ survey, 
with a median redshift of ~ 0.5 and a median mass of 
-^500 ~ 2.3 X 10^"^ M0//1. This sample is representative 
of the cluster yield and properties for the full survey, 
which was completed in November 2011 and has detected 
~500 clusters. Therefore it is useful to consider how the 
method used in this work will be applied to the full sur- 
vey, and what level of improvement we can expect on the 
cosmological constraints. 

Using the SPTcL+i?o+BBN data set, we found that 
our constraints are currently limited by both statistical 
uncertainty and the SZ-Yx scaling uncertainty. Both 
uncertainties would be improved by adding more SPT 
clusters with additional Yx measurements. Recently the 
SPT collaboration was awarded a Chandra X-ray Vision- 
ary Project (XVP) to complete X-ray observations of the 
80 most significant clusters at z > 0.4 detected in the first 
2000 deg^ of the SPT-SZ survey. As argued in Section 
5.1.2, we would need > 50 clusters with Yx measure- 
ments for the C — calibration to be limited by the 
Yx — M500 uncertainty. With this many clusters, the 
statistical uncertainty on w should decrease to a level 
below the systematic uncertainty from the X-ray scal- 
ing relation, Sw = ±0.15. Combining the full 2500 deg^ 
SPT-SZ survey with the Chandra XVP observations, we 
would be limited to this constraint from the current cal- 
ibration of the Yx — M500 relation. 

To reduce the systematic uncertainty further, we would 
need more accurate cluster mass estimates than currently 
exist from X-ray measurements alone. In Section 6.2, 
we found that the X-ray scaling systematics were cur- 
rently limited by the uncertainty in Ax and Cx, whose 
fractional uncertainty was 10% and 50%, respectively. 
Reducing their uncertainty by a factor of two, would re- 
duce their contribution to the systematic uncertainty to 
Sw = ±0.037 and ±0.074, respectively. This would effec- 
tively correspond to an overall mass calibration uncer- 
tainty of 5% with an additional 6% uncertainty in the 
evolution of the mass calibration between z = 0.0— 1.1. 

This level of mass calibration should be achievable by 
incorporating additional data sets, such as optical veloc- 
ity dispersion (White et al. 2010) or weak lensing mea- 
surements (Hoekstra 2007; Becker & Kravtsov 2011). For 
example, in massive clusters, the scatter in weak lens- 
ing mass estimates is expected to be ~20% (Becker & 



Kravtsov 2011). Therefore, with weak lensing observa- 
tions of two sets of ^15-20 clusters at low and high- 
redshift, this level of accuracy should be achievable. To- 
wards this goal, the SPT collaboration has been approved 
for weak lensing observations of ^^35 SPT-detected clus- 
ters spanning 0.30 < z < 1.3 using the Magellan and 
Hubble telescopes. Additionally, the SPT collaboration 
has been approved for optical velocity dispersion observa- 
tions of ~100 SPT-detected clusters using the Very Large 
Telescope (VLT) and a large NOAO program on Gem- 
ini South. With these data sets, we expect to achieve 
the factor of two improvement in mass calibration, as 
discussed above. 

Applying this calibration to the full 2500 deg^ 
SPTcL+-f^o+BBN data set, we should constrain w with 
an accuracy of ~8%, or a factor of ~4.5 tighter than 
the current SPTcl+^^^o+BBN constraints. This im- 
proved constraint would be comparable to the current 
constraints from the CMB+BAO+SNe data, and would 
be an independent systematic test of the standard dark 
energy paradigm by measuring the effect of dark en- 
ergy on the growth of structure. Combining the existing 
CMB+BAO+SNe data with the 2500 deg^ SPT clus- 
ter sample, the uncertainty from the SZ-Yx scaling and 
the chister sample size is expected to be negligible com- 
pared to the uncertainty contributed by the improved 
cluster mass calibration. In this case, the SPT cluster 
data would contribute an uncertainty of only ^ 1% to 
the significantly improved constraint on w from the com- 
bined data set. 

8. CONCLUSIONS 

We use measurements from the SPT-SZ cluster survey 
in combination with X-ray measurements to constrain 
cosmological parameters. We have described and imple- 
mented a method that simultaneously fits for cosmolog- 
ical parameters and the scaling of the SZ and X-ray ob- 
servables with cluster mass. The method is generalizable 
to multiple cluster observables, and self-consistently ac- 
counts for the effects of cluster selection and uncertainties 
in cluster mass calibration on the derived cosmological 
constraints. We apply this method to a SZ-selected cat- 
alog of 18 galaxy clusters identified in 178 deg^ of the 
2500 deg2 SPT-SZ survey. This is the first analysis of 
an SZ survey to directly incorporate X-ray observations, 
which has reduced the uncertainty on both the cluster 
mass calibration and the cosmological constraints. 

For a ACDM cosmology, we find that 
the SPTcL+i?o+BBN data best constrain 
(78(1^™/0.25)°-30 = 0.785 ± 0.037, where the total 
uncertainty consists of an approximately equal amount 
of statistical and systematic uncertainty. These con- 
straints are consistent, and comparable, with other 
constraints using X-ray-selected (Vikhlinin et al. 2009b; 
Mantz et al. 2010c) and optically-selected (Rozo et al. 
2010) cluster samples. In combination with measure- 
ments of the CMB power spectrum from the SPT data 
and the seven-year WMAP data, the SPT cluster data 
constrain as = 0.795 ± 0.016 and f2„ = 0.255 ± 0.016, a 
factor of 1.5 improvement on each parameter over the 
constraints from the CMB data alone. 

We consider several extensions beyond a ACDM cos- 
mological model by including the following as free param- 
eters: the dark energy equation of state {w), the sum of 
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the neutrino masses (Errii/), the effective number of rela- 
tivistic species (Nes), and a primordial non-Gaussianity 
(/jvl). 

For a wGDM cosmology, the SPTcL+-ffo+BBN data 
constrain w = -1.09 ± 0.36 and ag = 0.773 ± 0.088, 
consistent with dark energy being due to a cosmo- 

logical constant, and with comparable uncertainties to 
constraints from the CMB data alone. Using the 
CMB+BAO+SNc+SPTcL data set, we constrain w = 
-0.973 ± 0.063 and as = 0.793 ± 0.028, a factor of 1.25 
and 1.4 improvement, respectively, over the constraints 
without SPT cluster data. The uncertainty on w consists 
of approximately equal contributions from statistical un- 
certainty, systematic uncertainty from SNe, and system- 
atic uncertainty from cluster scaling relations, with the 
latter contributing an uncertainty of 5w = ±0.023. 

We next consider a ACDM cosmology with a non- 
zero neutrino mass. Using a CMB+iJo+BAO+SPTcL 
data set, we constrain the sum of the neutrino masses 
Em^ to be < 0.33 eV at 95% confidence, a factor of 1.4 
improvement over the constraints without SPT cluster 
data. We find even tighter constraints when we exclude 
the BAO data set, which tend to favor a lower value 
of Hq and therefore a higher neutrino mass. Using a 
CMB+i?o+SPTcL data set, we constrain 'Em^< 0.28 eV 
at 95% confidence. We also consider a model with a free 
effective number of relativistic species, A^effi to explain 
the increased damping that is observed in the the CMB 
power spectrum. Using a CMB-f-i?o+BAO+SPTcL data 
set, we jointly measure Ngg= 3.91 ± 0.42 and J^m^= 
0.34±0.i7 eV, while constraining Sm^< 0.63 eV at 95% 
confidence. 

Finally, we consider a ACDM cosmology where we al- 
low the number of observed clusters to be affected by 
non-Gaussian density fluctuations characterized by the 
parameter Jnl- Using a CMB+SPT^l data set, we 
measure fNL= —192 ± 310, consistent with zero non- 
Gaussianity. 

The results presented in this paper use 18 clusters 
from 178 deg^ of the 2500 deg^ SPT-SZ survey, and are 
limited by the combination of the cluster sample size 
and mass calibration. The SPT-SZ survey was com- 
pleted in November 2011, and has detected ^500 clusters 
with a median redshift of ^ 0.5 and a median mass of 
Afsoo ~ 2.3 X IO-^^Mq/Zi. Ongoing X-ray, weak lensing, 
and optical velocity dispersion observations of SPT-SZ- 
selected clusters will be used to produce an improved 
cluster mass calibration of the sample. The full SPT-SZ 
survey and improved mass calibration will produce con- 
straints on w comparable to current constraints from the 
combination of CMB-hBAO-f-SNe data, and would rep- 
resent an independent systematic test of the standard 
dark energy paradigm by measuring the effect of dark 
energy on the growth of structure. The combination of 
CMB-hBAO-hSNe data with the SPT cluster sample will 
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APPENDIX 

X-RAY OBSERVATIONS AND RESULTS 

In Table 7, we give updated X-ray observables for the clusters used in this work, as discussed in Section 2.1.3. For 
the five clusters without new measurements, we give the results from All directly, in order to provide a complete 
listing for the cluster sample. In Table 8, we give the complete list of Chandra observation identifications (ObsIDs) 
used for clusters with new Chandra observations 
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Note. — X-ray observables for clusters with Chandra or XMM-Newton observations. For clusters with new spectroscopic 
redshifts or new Chandra obser-vations, we have recalculated their X-ray observables, as described in Section 2.1.3. To maintain 
consistency with All, all X-ray observables are calculated assuming a preferred ACDM cosmology using WMAP7+BAO±i?o 
data with Qm = 0.272, = 0.728 and Ho = 70.2 km s"' Mpc"' (Komatsu et al. 2011). 

''X-ray observables taken from All. 
''Updated for new Chandra observations. 
■^Updated for new spectroscopic redshift. 



TABLE 8 

Clusters with new Chandra X-ray Observations 





All 


This Work 




Name 


z 


Exposure 


Source 


Exposure 


Source 


ObsIDs 






[ks] 


Counts 


[ks] 


Counts 




SPT-CL J0528-5300 


0.765 


36.5 


356 


115.9 


1732 


9341, 10862, 11996, 11747, 11874, 12092, 13126 


SPT-CL J0533-5005 


0.881 


41.5 


201 


67.7 


344 


11748, 12001, 12002 


SPT-CL J0546-5345 


1.067 


55.6 


1304 


67.8 


1512 


9332, 9336, 10851, 10864, 11739 


SPT-CL J0551-5709 


0.423 


19.8 


876 


33.2 


1212 


11871, 11743 


SPT-CL J2359-5009 


0.774 


57.9 


713 


122.4 


1522 


9334, 11742, 11864, 11997 



Note. The ObsIDs refer to all the Chandra observations used in this work. The ObsIDs that arc new, relative to All, are 
higiilighted in bold. For clusters not listed here, we use the same Chandra and XMM-iVewton observations as listed by All. 



LIKELIHOOD MODIFICATION TO ACCOUNT FOR COSMOLOGICAL DEPENDENCE OF Yx 

In Section 3.2 wc presented a procedure for translating the theoretical mass function, dN /dMdz, into observable 
space, dN/d^dYxdz. Under the assumption that this transformation is independent of the cosmological and scaling 
relation parameters p we are ultimately trying to recover, this procedure modifies the log likelihood by a constant 
offset. 

However, in the case of Yx, this assumption ceases to hold true as Yx is a derived quantity, calculated explicitly for 
each new value of p (i.e., at each likelihood evaluation in the MCMC). This added subtlety can be addressed in the 
following fashion: Let us define Y^ as the value of Yx when evaluated at some reference point in parameter space p* . 
In order for the probability contained in a differential volume to be independent of a change of variables, we need to 
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multiply by the Jacobian of the transformation, as follows 

where is the ratio of the Yx value of the i*^ cluster at p to its value at p*. Expressing this in terms of log 

probability, and ignoring constant offsets, we obtain 

In P{z, I \p) = In P{z, I Y: Ip) + ^ In F,,^ (B2) 

i 

This results in the straightforward prescription of adding ^InYx to the likelihood at each step in the MCMC, a 
process very similar to that suggested and performed in Mantz et al. (2008) and Vikhlinin et al. (2009a). 

MASS ESTIMATES 

We present posterior mass estimates for all 18 clusters considered in this work in Table 9. Where applicable, these 
are joint X-ray and SZ posterior mass estimates, for clusters without X-ray data we use the SZ posterior mass estimate. 
We calculate a probability density function on a mass grid at each point in the ACDM chain that was calculated using 
the CMB+SPTcL data, from Section 4.1. The probability density functions arc combined to obtain a mass estimate 
that has been fully marginalized over all cosmological and scaling relation parameters. We report the mean and the 
68% confidence interval for the mass estimate. 

TABLE 9 

Mass Estimates for the SPT Cluster Catalog 



Object Name (, z M6oo(PcHt)(10" Mq/i^q^) 



SPT- 


CL 


J0509-5342 


6.61 


0.463 


5 


11 


± 





68 


SPT- 


CL 


J0511-5154'' 


5.63 


0.74 


3 


36 


± 





86 


SPT- 


CL 


J0521-5104'' 


5.45 


0.72 


3 


21 


± 





86 


SPT- 


CL 


J0528-5259 


5.45 


0.765 


2 


96 


± 





54 


SPT- 


CL 


J0533-5005 


5.59 


0.881 


2 


54 


± 





54 


SPT- 


CL 


J0539-5744'' 


5.12 


0.77 


2 


93 


± 





86 


SPT- 


CL 


J0546-5345 


7.69 


1.067 


4 


79 


± 





64 


SPT- 


CL 


J0551-5709 


6.13 


0.423 


3 


61 


± 





54 


SPT- 


CL 


J0559-5249 


9.28 


0.611 


6 


36 


± 





79 


SPT- 


CL 


J2301-5546'' 


5.19 


0.748 


3 


00 


± 





86 


SPT- 


CL 


J2331-5051 


8.04 


0.572 


4 


89 


± 





68 


SPT- 


CL 


J2332-5358 


7.30 


0.403 


6 


21 


± 





79 


SPT- 


CL 


J2337-5942 


14.94 


0.781 


7 


68 


± 


1 


04 


SPT- 


CL 


J2341-5119 


9.65 


0.998 


5 


14 


± 





71 


SPT- 


CL 


J2342-5411 


6.18 


1.074 


2 


75 


± 





46 


SPT- 


CL 


J2355-5056 


5.89 


0.320 


3 


96 


± 





54 


SPT- 


CL 


J2359-5009 


6.35 


0.774 


3 


32 


± 





54 


SPT- 


CL 


JOOOO-5748 


5.48 


0.701 


4 


04 


± 





68 



These clusters have only SZ data, and no X-ray observations. 
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